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Molecular Phylogeny and Taxonomic Revision of the Woolly Lemurs, 

Genus A vahi (Primates: Lemuriformes) 


RambinintsoaAndriantompohavana, RunhuaLei, John R. Zaonarivelo, Shannon E. Engberg, 
Gerard Nalanirina, Susie M. McGuire, Gary D. Shore, Justin Andrian.asolo, 

Kelly Herrington, Rick A. Brenneman, and Edward E. Louis, Jr. 

Abstract 

Approximately 3,000 base pairs of mitochondrial DNA sequence data (control region or 
D-loop, along with a fragment including a partial segment of cytochrome oxidase subunit III gene 
to Leucine-tRNA gene), morphological and phenotypic characters were used to investigate the 
phylogenetic relationships among the recognized woolly lemurs (genus Avahi ) and to other lemur 
genera of Madagascar. The molecular data support the previously recognized Avahi species, but 
it has also revealed additional unique biodiversity. Here, we present further taxonomic revisions 
of the genus Avahi, by elevating the two recently described subspecies, A. m. meridionalis and 
A. m. ramanantsoavana , to species level, along with a description of one new species. We also 
conducted the first phylogenetic analysis of all recognized species of woolly lemurs, all repre¬ 
sentatives of the family Indriidae, and present the relationships between the recognized species 
of the genus Avahi and this newly described species. These results underscore the urgency to 
initiate further detailed studies in previously unstudied sites throughout Madagascar in order to 
better define lemur species. 

Key words: Avahi , D-loop, Madagascar, prosimian, systematic, woolly lemurs 


Introduction 


Due to its unique species biodiversity and the con¬ 
tinued pressure from human encroachment in the form 
of habitat destruction and fauna and flora utilization, 
Madagascar has been placed at the top of conservation 
priority lists, or hotspots (Myers 2000). Distributed 
throughout the island, the prosimians of Madagascar 
are particularly susceptible to extinction risks due to 
their relatively small, fragmented geographic ranges, 
and the potential loss through edge effects and constant 
habitat loss (Jernvall and Wright 1998; Reed 2004; 
Ezard and Travis 2006). Since 2000, multiple molecu¬ 
lar genetic and morphological studies within the genera 
Avahi , Cheirogaleus, Lepilemur, Microcebus, and Mina 
have led to a significant increase in the number of 
recognized species (Zimmermann et al. 1998; Groves 
2000; Rasoloarison et al. 2000; Kappeler et al. 2005; 
Thalmann and Geissmann 2005; Andriaholimrina et 
al. 2006; Andriantompohavana et al. 2006; Louis et 
al. 2006a, 2006b; Zaramody et al. 2006). With the 


addition of newly described species, the once broad 
ranges of the previously recognized nocturnal species 
have dramatically decreased and concurrently made 
the risk of extinction of both recognized and newly 
described lemurs much greater (a smaller distribution 
and a diminished number of represented subpopula¬ 
tions). Currently, all lemurs are protected under the 
Convention oflntemational Trade of Endangered Spe¬ 
cies (CITES) and are designated by thelUCN/SSC Red 
List Categories from critically endangered to threatened 
(IUCN 2004), Recent revisions underscore the need to 
resolve the taxonomy and phylogeny of lemurs so that 
a scientifically rational approach to their conservation 
and management can be developed and implemented. 
Furthermore, the numerous taxonomic revisions of 
lemurs have led to subsequent modifications of histori¬ 
cal distributions and population estimates that require 
a persistent re-evaluation of conservation protection 
status (Martin 2000; Louis et al. 2006b). 
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The woolly lemurs (genus Avahi), a group of 
strictly nocturnal primates that are regarded as pheno- 
typically indistinct with varying degrees of brown, red, 
white, and gray pelage, are a prime example of species 
whose taxonomic/evolutionary diversity until now has 
been largely unrecognized and underrepresented. Tra¬ 
ditionally, two species of Avahi have been recognized, 
one from eastern Madagascar, Avahi laniger , and one 
from western Madagascar, Avahi occidental is (Fig. 
1). Initially, both species were considered subspecies 
of Avahi lanigei\ but morphological, molecular, and 
cytogenetic differences between the eastern woolly 
lemur, Avahi laniger, and the western woolly lemur, 
Avahi occidentalism have conclusively established spe¬ 
cies level variance (Rumpler et al. 1990; Thalmann and 
Geissmann 2000; Zaramody et al. 2006; Figs. 2-4). 

The potential for additional species within the 
genus Avahi was predictable, given the discrete and 
often fragmented distributions, the scattered or spotty 
occurrence within these fragments, and the widespread 
area that the genus Avahi is found especially on the 
eastern coast of Madagascar (Martin 2000; Thalmann 
and Geissmann 2000). A comprehensive study to 
define these populations is critical above and beyond 
the threat of habitat loss because the highly foli- 
vorous dietary requirements of the woolly lemur can 
bring about a selective liability that can potentially 
limit habitat suitability. Consequently, Thalmann and 
Geissmann (2000. 2005) based on morphological and 
vocalization data described two additional woolly 
lemurs from western Madagascar, A. unicolor and A. 
cleesei , from the Ampasindava region and Tsingy de 
Bemaraha, respectively (Figs. 1, 5-7). Furthermore, 
Zaramody et al. (2006) described two new species 
from the southeastern coast, A. peyrierasi (south of 
the Mangoro/Onive Rivers to Ranomafana National 
Park) and A. meridionals (Manombo Special Reserve 


to Andohahela National Park; Fig. 1). Avahi meridi¬ 
onals was further subdivided into two subspecies, A. m. 
meridionals and A. m. ramanantsoovana, from Sainte 
Luce/Andohahela National Park and Manombo Special 
Reserve, respectively (Fig. 1). Here, we present further 
taxonomic revisions of the genus Avahi , along with a 
description of one new species. We also conducted 
the first phylogenetic analysis of all recognized spe¬ 
cies of woolly lemurs, all representatives of the family 
Indriidae, and present the relationships between the 
recognized species of the genus Avahi and this newly 
described species. 

As previously discussed in Andriantompohavana 
et al. (2006), Louis et al. (2006a, 2006b), and Thalmann 
and Geissmann (2005), the utilization of whole vouch¬ 
ers as the designated holotvpe for a new species is not 
a prerequisite, but that opportunistic collections can 
later supplement morphological, vocalization, and/or 
molecular data in combination with curated blood 
and/or tissue samples. The woolly lemurs are a prime 
candidate for this methodology because the highly foli- 
vorous dietary requirements of this group of lemurs cur¬ 
rently precludes any attempts to curate t; live vouchers” 
(Thalmann and Geissmann 2005). Total genomic DNA 
for the three paratype specimens are currently curated 
at the Museum of Texas Tech University (TK125757; 
TK125758; TK 125759). Additionally, an electronic 
database that includes all Avahi field data and photo¬ 
graphs, including data for the paratype specimens, is 
curated at the Museum of Texas Tech University. The 
database is stored in the Type Specimen Collection in 
multiple media formats. This collection of field data 
and photographs, as well as additional tables and fig¬ 
ures, also are available online at the website of Omaha’s 
Henry Doorly Zoo. See Appendix I for a directory of 
appropriate website addresses. 


(text continued on page 10) 



Andrianiompohavana et al.— Woolly Lemurs of Madagascar 



Mahavavy Nord- 

Antafondro- 

Sambirano- 

0 100 


Mahajamba 

Mariarano 

Bot&iboka 


Mahavavy $ud 
Maningoza 


Manambaho 
Manambolo 

Tsingy de 
Bemaniha 

Tsiribihina 

Morondava 


Mangdky 


FihcrLTiana 

Onilahy 


Lok) 


Maevarano 


Sofia 


Bemnrivo 


Antainambalana 

Mananara-Nord 


Maningoiy 
— On She 


Mantadia 
Maromizaha 

Mango i;o 

Fandriana 

Namorurta 

Ranomafana 


Lima 


/ 

Monaraiidra 


Paraony 

Maiiamp&tnma 

Manombo 

Mananara 

Andohahcla 

(Manangotry) 


Mandril re 


0* Aw hi peyriefttsi 

nr. taman&nisnavums 

0 Avtiftim- wcrirffonatis 
Availi lonigvr 
Avaki anicoittr 
Aviihl cicesci 
} Avail f ucvidvii ta fin 

< Foul eoSor of site denotes 
we! coast Avahi spec its. 
alt formerly Avtthi occitlcniulis 
^ FimI color of situ denotes 
cast coast Avaki species, 
all formerly Avaki tamper 


Figure 1. Sample distribution map of the woolly lemurs of Madagascar. 
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Figure 2. Avahi occidentalis , western woolly lemur, at Mariarano Classified Forest. Photo by Edward E. 
Louis, Jr. 
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Figure 3. Avahi occidentalism western woolly lemur, at Mariarano Classified Forest. Photo by Edward E. 
Louis, Jr. 
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Figure 4. Avahi laniger, western woolly lemur, at Mananara-Nord National Park (Ivontaka-Sud). Photo by 
Edward E. Louis, Jr. 
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Figure 5. Avahi unicolor, Sambirano woolly lemur, at Antofondro Classified Forest (Maromiandra). Photo by 
Edward E. Louis, Jr. 
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Figure 6. Avahi unicolor , Sambirano woolly lemur, at Antofondro Classified Forest (Maromiandra). Photo 
by Rambinintsoa Andriantompohavana. 





Andrianiompohavana et al.— Woolly Lemurs of Madagascar 


9 



Figure 7. Avahi cleesei, Cleese’s woolly lemur, at Tsingy de Bemaraha National Park. Reprinted with permission of the 
author and journal (Thalmann and Geissmann 2006). Photo by Urs Thalmann. 
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Methods 


Sampling .—All lemurs investigated in this study 
were wild-caught (sampling localities are depicted in 
Fig. 1) and were immobilized with a CO, projection 
rifle or blowgun with 10 mg/kg of Telazol (Fort Dodge 
Animal Health; Overland Park, Kansas). Four 2.0 mm 
biopsies and 1.0 cc per kilogram of whole blood were 
collected from each sedated animal and immediately 
stored in a room temperature storage buffer (Longmire 
etal. 1992). The location of the immobilized lemur 
was recorded using a global positioning system (Ap¬ 
pendix la). Genomic DNA was extracted from a 2.0 
mm ear punch using a phenol-chloroform extraction 
(Sambrook et al. 1989; Appendix II). Beginning in 
2000, while the animal was sedated, we placed a 
HomeAgain microchip (Schering-Plough Veterinary 
Corp.; Kenilworth, New Jersey) subcutaneously be¬ 
tween the scapulae of each lemur (Appendix II). This 
procedure was used to field catalog each animal with 
a unique recognition code to provide for the capabil¬ 
ity to positively re-identify all captured individuals 
during any future immobilizations. Measurements 
were taken on the sedated animals as follows, weight 
(+/- 0.01 kilograms); head crown ((total length from 
tip of the nose (not the soft tissue) to the occipital 
crown) +/-0.1 cm); body length ((total length of body 
from the occipital crown of the head to the base of 
tail) +/-0.1 cm); tail length ((total length from base of 
tail to the end of the last caudal vertebra) +/-(). 1 cm); 
forelimb measurements (thumb (total length of the 
thumb (proximal, middle, and distal phalanges) from 
the distal tip to proximal point of the last bone of the 
thumb (claw not included) +/-0.1 cm), longest digit 
(total length of the phalange only (proximal, middle, 
and distal phalanx, claw and metacarpal not included) 
+/-0.1 cm), and hand (total length of hand (phalanges 
and metacarpals included to the carpal joint) +/-0.1 
cm)); radius/ulna ((carpal joint (styloid process of the 
radius) to the olecarnon tip of the ulna) +/-0.1 cm); 
humerus ((greater tubercle of humerus to the end of the 
lateral condyle of the humerus) +/-0,1 cm); hindlimb 
measurements ((toe (total length of the toe (proximal, 
middle, and distal phalanges not including claw) from 
tip to insertion point of the last bone of the hind toe 
(1st digit) +/-0.1 cm), longest digit (total length of 
the phalange only (claw and metatarsal not included) 
+/-0.1 cm), foot (total length of hand (phalanges and 
metacarpals included to the tarsal joint, not including 


claw) +/-0.1 cm), tibia ((calcaneal tuber to the proxi¬ 
mal tibial tuberosity ) +/-0.1 cm), and femur ((greater 
trochanter of the femur to the distal point of the lateral 
condyle of the femur) +/-0.1 cm)); testes ((width and 
length of the right and left testes in males) +/-0.1 mm); 
and upper and lower canines ((length of the upper and 
lower right canines from the tip to the gum line) +/-0.1 
mm). For presentation purposes, we present average 
measurements and standard deviations of the weight, 
head crown, body length, tail length, forelimb, and 
hindlimb in this publication following the guidelines 
of Smith and Jungers (1997; Table 1A-B). All other 
measurements, e-voucher photographs, and data are 
available at http://www.omahazoo.com/ccr/index. 
asp?page=/ccr/genetics/genhome.htm (Appendix la). 
All traits were tested for significance among species 
using Analysis of Variance (ANOVA) performed using 
an online statistical service (Kirkman 1996). Species 
means were plotted in histogram form with the appro¬ 
priate 95% confidence intervals (Appendix Id). 

Data Collection .—To correlate our data with 
previously published molecular studies, we analyzed 
the following regions of the mitochondrial DNA 
(mtDNA): the displacement loop or control region 
(D-loop; Baker et al. 1993; Wyner et al. 1999) and a 
fragment of the cytochrome oxidase subunit III gene 
(COIII), NADH-dehydrogenase subunits 3, 4L, and 
4 (ND3, ND4L, and ND4), as well as the tRNA Gly , 
tRNA Arg , tRNA His , tRNA Ser , and partial tRNA Leu genes 
(subsequently referred to as the PAST fragment; 
Pastorini et al, 2000, 2003). Using 50 nanograms of 
genomic DNA, the D-loop (562-563 base pairs (bp)) 
and the PAST (2380 bp) fragments were amplified 
using the following conditions: 94° C for 30 sec, 47° 
C for 45 sec, 72° C for 45 sec for 34 cycles. Since 
potential nuclear insertions or mitochondrial pseu¬ 
dogenes within the nuclear genome can be amplified 
inadvertently, we chose to minimize this likelihood 
by amplifying both mitochondrial DNA regions as 
intersecting or overlapping segments (Zhang and 
Hewitt 1996). Consequently, the PAST fragment was 
generated from six amplified segments. Additionally, 
to further eliminate amplification of nuclear inser¬ 
tions, a technique that is species independent and 
both rapid and efficient derived from the degenerate 
oligonucleotide-primed PCR method (DOP-PCR; 
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Telenius et al. 1992) was used to generate the PCR 
products. Adapting this LL-DOP-PCR (long products 
from low quantity ), the sequence data were generated 
for the D-loop fragments and PAST sequence gener¬ 
ated from overlapping segments were confirmed. The 
samples were electrophoresed on a 1.2% agarose gel to 
verify the PCR product and purified using QIAquick 
PCR purification kit (Qiagen; Valencia, California). 
Using the BigDye terminator cycle sequencing ready 
reaction kit by Applied Biosystems the sequence 
was analyzed by capillary electrophoresis with an 
Applied Biosystems Prizm 3100 genetic analyzer 
(Applied Biosystems, Inc.; Foster City, California). 
A suite of internal sequencing primers from Pastorini 
et al. (2000, 2001) were used to generate the PAST 
fragment. Additionally, PCR and sequencing prim¬ 
ers specific for genus Avahi were designed for the 
D-Loop fragment, and PCR and sequencing primers 
specific for ND4 region were utilized for the PAST 
fragments (Arevalo et al. 1994; Appendix lb). The 
sequence fragments were aligned to generate a con¬ 
sensus sequence using Sequencher (Gene Corp; Ann 
Arbor, Michigan), and the consensus sequences were 
aligned using ClustalX (Thompson et al. 1997). All 
aligned sequences are available from the first author 
upon request. All sequences have been deposited 
in GenBank and the sequence data and information 
are available from the referenced accession numbers 
(Appendix 11). 

Phylogenetic Analysis. —Maximum-parsimony 
analyses (MP) were performed for the phylogenetic 
study of the D-loop, PAST, and combined (D-Loop 
and PAST fragments) sequence data with PAUP Ver¬ 
sion 4.0b 10 software (Swofford 2001). Heuristic 
searches were completed using the random addition 
sequence (1000 replicates) with the tree bisection- 
reconnection branch swapping routine. Gaps were 
considered as a fifth character in MP analyses, but 
were treated as missing data in the NJ analyses. Boot¬ 
strap analyses were accomplished with 1000, 3000. 
and 4000 pseudoreplicates with the D-loop, PAST, 
and combined data sets, respectively. Only nodes 
with greater than 50% support were reported. The 
maximum likelihood (ML) analyses were performed 
using the program PAUP Version 4.0b 10 software 
(Swofford 2001). Due to the large number of taxa and 


characters and the resulting computational intensity, 
we pruned the combined sequence dataset by choos¬ 
ing taxa representing distinct haplotypes supported 
in an initial neighbor-joining (NJ; Saitou and Nei 
1987) analysis. 

Bayesians inference analyses were conducted 
using MrBayes 3.0b4 (Huelsenbeck and Ronquist 
2001; Ronquist and Huelsenbeck 2003). The model 
of evolution was selected by using MrModeltest 2.2, a 
modified version of Modeltest 3.6 (Posada and Cran¬ 
dall 1998; Nylander 2004). A Markov Chain Monte 
Carlo (MCMC) run with four simultaneous chains 
and 1,000,000 generations was performed. Every 
hundredth generation, the tree with the best likelihood 
score was saved, resulting in 4000 trees. The 4000 
trees were condensed in a majority rule consensus 
tree using PAUP Version 4.0b 10 software (Swofford 
2001). Branch supports were assigned as posterior 
probabilities on the consensus tree. The pattern of 
sequence evolution was estimated by conducting a 
minimum spanning network generated with the pro¬ 
gram NETWORK Version 4.11 (Gonzales et al. 1998; 
Bandelt et al. 1999; Forster etal. 2001) and Arlequin, 
Version 2.0 (Schneider et al. 2000). 

In addition to character-based phylogenetic 
analysis of DNA sequences, PAUP software (Swof¬ 
ford 2001) was also used to calculate uncorrected 
pairwise distances for D-Loop, and PAST combined 
fragments ( c p ? ). As described in Davis and Nixon 
(1992), Wyner et al. (1999), Mayor et aL (2004), 
Andriantompohavana et al. (2006), and Louis et al. 
(2006a, 2006b) we utilized MacClade 3.01 (Maddison 
and Maddison 1992) and MEGA Version 2.0 (Kumar 
et al. 1993) in a diagnostic search to designate evo¬ 
lutionary significant units (ESU) using population 
aggregate analysis (PAA) of the D-loop (562-563 bp) 
and PAST (2380 bp) sequence data for Avahi. With 
the sequential addition of each individual without 
an a priori species designation, a PAA distinguishes 
attributes or apomorphic characters according to the 
smallest definable unit (Davis and Nixon 1992; Mayor 
et al. 2004; Andriantompohavana et al. 2006; Louis 
et al. 2006a, 2006b). 
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Results 


Following ANOVA among the species, sub¬ 
species and two additional types, only the tail length 
was found to be nonsignificant. The body length 
was slightly significant (P<0.05), the head crown, 
humerus and femur lengths were moderately signifi¬ 
cant (P<0.01), while all other measurements were 
strongly to highly significant (PO.OO 1). Means and 
standard deviations are reported in Tables 1A and 1B; 
histograms defining the 95% confidence intervals 
for all measurements are located in Appendix Id. 

Mitochondrial DNA sequence data were com¬ 
pleted for two fragments, D-loop and PAST fragment 
(approximately 2380 bp) for 98 individuals, repre¬ 
senting all six recognized species of woolly lemurs 
from a total of twelve sites (Fig. I; Appendix II). All 
new mtDNAsequences generated for this study were 
deposited in GenBank and can be acquired through 
the accession number (Appendix la: Appendix II). 
The sequence alignments for the data sets are avail¬ 
able from the first author upon request. The PAST 
fragment consists of the 3' end of the COI1I gene (30 
bp), the complete NADH-dehvdrogenase subunits 
ND3 (348 bp), ND4L (297 bp), and ND4 (1378 bp), 
along with the tRNAgenes. glycine (73 bp), arginine 
(73 bp), histidine (70 bp), serine (65 bp), and the 
5 5 portion of leucine (47 bp). The polyadenylation 
of COIII and ND4 genes, insertion of base pairs 
between ingroup/outgroup comparisons, and other 
alignment characteristics between lemurs and Homo 
are consistent with Pastorini et al. (2000). 

Based on the phylogenetic inferences oftheNJ, 
MP, and ML analyses of three sequence alignments 
(D-loop, PAST, and D-Loop-PAST fragment com¬ 
bined), two major Avahi subgroups are represented, 
differentiating the six recognized woolly lemur spe¬ 
cies (Figs. 8-14). The first subgroup corresponds 
to the western woolly lemurs, A occidentalism A. 
cleesei , and A. unicolor (Fig. 8). The second sub¬ 
group includes the eastern woolly lemurs, A. me¬ 
ridional is meridional is m A. m. ramanantsoavcina, A. 
peyrierasi , and A, laniger (Fig. 8). Two haplotype 
groups or types of A. peyrierasi were represented 


within the eastern woolly lemurs (Figs. 8-14). In 
addition, another group of woolly lemurs from the 
Fandriana region formed a distinct clade,zL species 
nova # 1, from the other recognized species (Fig. 8). 
There is high bootstrap support for the MP and NJ 
analysis with respect to the topology of the genera 
and species as every major node between the species 
clusters has a greater than 50% bootstrap support 
(Figs. 8, 11, and 13). The results from the popula¬ 
tion aggregate analysis of the D-loop and PAST se¬ 
quence data are presented in Tables 2A-B and 3A-B, 
respectively. Multiple diagnostic characters define 
each Avahi species, along with two types of Avahi 
laniger (only the results from the PAST fragment 
PAA designates attributes specific for each region), 
two types of Avahi peyrierasi , and the one newly 
described species, A. species nova #1 (Figs. 15-17; 
Tables 2 and 3). A review of the morphometric data 
for the six recognized species of woolly lemurs are 
summarized in Table 1 and are presented in Appendix 
Id, along with the two types of A, peyrierasi , and 
the newly described species (Appendix la; detailed 
morphological measurements of the individual 
animals are available at http://www.omahazoo. 
com/cer/index.asp?page=/ccr/genetics/genhome. 
htm) . The uncorrected *p' distance and the absolute 
genetic differences are presented in Table 4A-B. The 
highest average uncorrected pairwise distances and 
absolute genetic differences for the D-Loop sequence 
alignment were between the western and eastern 
woolly lemurs, ranging from 12.28 to 8.84 and 65 
to 47 (A. cleesei to A. peyrierasi type #2 and A. oc¬ 
cidental is to A . peyrierasi #3), respectively. The 
highest average uncorrected pairwise distances and 
absolute genetic differences for the PAST sequence 
alignment were between the western and eastern 
woolly lemurs, ranging from 11.72 to 10.20 and 251 
to 220 (A. occidental is to A. in. meridionalis and A. 
occidentalis to A . peyrierasi #3), respectively. The 
highest and lowest average uncorrected pairwise 
distances and absolute genetic differences for the 
D-Loop sequence alignment between the eastern 
woolly lemurs range from 5.99 to 3.43 and 43 to 
19 (A. m. meridionalis to A. laniger and A. species 


(text continued on page 35) 
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Figure 8. Neighbor-joining phylogram derived from the D-loop DNA sequence data from the 
77 Avahi individuals with 26 out-group taxa. Species designated according to the distribution 
in the current literature (Mittermeier et al. 2006; *Thalmann and Geissmann 2000; *Thalmann 
and Geissmann 2005; **Zaramody et al. 2006). Values above branches indicate number of 
changes between nodes. Values within circles indicate support of bootstrap pseudoreplicates. 
Avahi catalogue numbers that are listed in Appendix II, but are not shown in this figure, are 
presented as haplotype information in Appendix lc. 
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Figure 9. Maximum parsimony phylogram derived from D-loop fragment sequence data from 77 Avahi individuals (one of 144 
most parsimonious trees). Values above branches indicate number of changes between nodes. Length = 1,271; Cl = 0.5382; 
RI = 0.8786; RC = 0.4728; HI = 0.4618. 
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Figure 10. Maximum-likelihood phylogram derived from D-loop sequence data from 29 Avahi haplotypes. 
The phylogram presented with branch lengths proportional to the number of changes (values specified 
on the branches). We obtained the maximum likelihood phylogram (-In likelihood = 6398.15) from the 
D-loop alignment from a transition/transversions ratio of 2.28 (k = 4.66) and y shape parameter of 0.54. 
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Figure 11. Maximum parsimony phylogram derived from PAST fragment sequence data from 77 Avahi individu¬ 
als (one of 53 most parsimonious trees). Values above branches indicate number of changes between nodes. 
Length = 3,551; Cl = 0.4816; RI = 0.8721; RC = 0.4200; HI = 0.5184. 
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Figure 12. Maximum-likelihood phylogram derived from PAST sequence data from 77 Avahi individuals. The 
phylogram presented with branch lengths proportional to the number of changes (values specified on the branches). 
We obtained the maximum likelihood phylogram (-In likelihood = 18,123.61) from the PAST alignment from 
a transition/transversions ratio of 5.94 (k = 12.31) and y shape parameter of 0.34. 
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Figure 13. Neighbor-joining phylogram derived from the D-loop and PAST combined DNA sequence data from the 39 Avahi 
haplotypes with 19 out-group taxa. Values above branches indicate number of changes between nodes. Values within circles 
indicate support of bootstrap pseudoreplicates. 
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Figure 14. Maximum parsimony phylogram derived from the D-loop and PAST combined DNA sequence data from the 39 
Avahi haplotypes with 19 out-group taxa (one of 20 most parsimonious trees). Values above branches indicate number of 
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Table 2 A. Diagnostic nucleotide sites from the D-Loop Pairwise Aggregate Analysis 
(PAA) of Avahi. _ _ 
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Table 2A. (cont.) 
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Table 3 A. Summary of Population Aggregate Analysis (PA A) D-Loop diagnostic sites for Avahi species. Refer to Table2A. *No 
character o r att ribute^ is avai lable for this fragme nt (A. betsi leo was formerly referred to as Avahi species nova Hi). 


Species 

Fragment Size (bp) 

PAA base pair location 

Avahi laniger type U 1 (Mantadia) 

562 

543 

Avahi laniger type #2 (Mananara-Nord) 

562 

* 

Avahi peyrierasi type #2 

562 

336, 548 

Avahi peyrierasi type #3 

562 

172. 164, 195,333, 502,505 

A vahi meridional is ramanat ttsoa vana 

562 

15S. 194. 199,525 

Avahi meridional is meridionals 

562 

16, 64. 159. 422,506 

Avahi betsileo 

562 

20, 227, 277, 309, 334, 348.521 

Avahi unicolor 

563 

106, 288, 313, 396.397.401 

Avahi cleesei 

563 

64. 112, 275, 276. 279,287, 300, 302. 332 

A.vahi o cc idea talis 

563 

295, 296, 314, 499 


Table 311. Summary of Popul ation Aggregate Analysis (PA A) PAST diagnostic si tes for Avahi species. Refer to Table! B . 

Fragment 

Species PA A base pair location. 

Size (bp) 


Avahi hunger type #1 (Mantadia) 2364 

Avahi laniger type #2 (Mananara-Nord) 2364 


Avahipeyrierasi type #2 2364 

Avahipeyrierasi type #3 2364 

Avahi meridionalis rantanantsoavana 2364 

Avahi meridionals meridiomdis 2364 

Avahi betsileo 2364 

Avahi unicolor 2364 

A vahi cle esei 2364 

Avahi occidentals 2364 


671,957. 1117, 1708 
576.835,1390,1397. 1708, 1870 

54, 108, 246. 701,962. 1408. 1584, 1861, 1946, 2096.2153, 1052 

14, 75. 155,344,572, 590,698,945, 1124, 1139, 1341. 1353, 1525, 1537, 1818, 1835, 1886, 
1960.2059,2111,2123 

29, 279, 938, 940, 972, 985, 1261, 1273, 1381, 1384, 1407, 1807 

85,289,291, 339,708, 847, 935, 952. 958. 1060. 1 I 14, 1131, 1270, 1336, 1339, 1612,1771, 
2013, 2355 

108.336, 838,920, 1052, 1053. 1363, 1813, 1837 

6,93, 114, 133,501,636.916, 1052, 1135, 1279, 1454, 1849, 1872,2032,2116,2131 
8, 118,213,258,522, 704,785,928, 1057, 1213, 1273, 1552, 1585, 1786 

47, 111, 121,500, 540, 665. 764, 783, 826, 1075, 1138. 1168, 1333, 1477, 1519, 1525, 2053, 
2059,2116,2138,2331 
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Figure 15. Avahi laniger. eastern woolly lemur, at Mantadia National Park. Photo by Edward E. Louis, Jr. 
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Figure 16. Avahipeyrierasi Type #2, Peyrieras’ woolly lemur (Type #2), at Ranomafana National Park (Am- 
batolahy Dimy). Photo by Edward E. Louis, Jr. 
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Figure 17. Avahipeyrierasi Type #3, Peyrieras’ woolly lemur (Type #3), at Ranomafana National Park (Ambatolahy 
Dimy). Photo by Edward E. Louis, Jr. 
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Table 4A. Genetic distance matrix for D-Loop sequence data for Avahi species. The numbers represent the following Avahi species: 
[}] Avahi peyrierasi type i'2; [2] Avahi peyrierasi type H3; [3] Avahi ramanantsoavana; [4] Avahi betsileo; [5] Avahi meridionali.s; 
[6] Avahi laniger: [7] Avahi unicolor; fSJ Avahi cleesei; and [9] Avahi occidental is. Genetic distance bused on absolute differences 
is displayed above the diagonal, and genetic distance based as a percentage is displayed below the diagonal (A, betsileo was formerly 



1 

2 

a 

D 

4 

5 

6 

7 

8 

9 

1 


3t 

29 

28 

00 

42 

62 

65 

59 

2 

5,28+0.95 


21 

19 

29 

43 

50 

56 

47 

j 

4.7S±0.87 

3.67*0.80 


24 

21 

42 

56 

59 

51 

4 

4,75+0.91 

3.43*0.77 

4,21 ±0,85 


28 

42 

58 

62 

53 

5 

5.56=0.98 

5.18*0.97 

3.51+0,75 

5.04+0,96 


43 

56 

58 

51 

6 

5.31 ±0.86 

5.64*0.93 

5.85±0.92 

5.50±0,9I 

5.99*0,97 


58 

63 

56 

7 

11.63*1 .43 

9.4 7± 1.26 

10.67*1.30 

11.30*1.43 

10.83*1,42 

9.61*1,23 


24 

15 

8 

12.28*1.46 

10.52*1.39 

11.22*1.42 

12.07*1.52 

11.13*1 43 

10.78*1.36 

4.17*0.85 


22 

9 

11.03±1.36 

8.84±1.22 

9.62+1.28 

I0.24±1.35 

9.75*1-34 

9.04*1.15 

2.76*0.70 

3.86*0,80 


Table 4B. Genetic distance matrix for PAST sequence data for Avahi species. The numbers represent the following Avahi species: 

[!] Avahi peyrierasi type H2; [2] Avahi peyrierasi type 43; [3] Avahi ramanantsoavana; [4] Avahi betsileo; [5] Avahi meridional is; 
[6] Avahi laniger; [7] Avahi unicolor; [8] Avahi cleesei; and [9] Avahi occidcnlalis. Genetic distance based on absolute differences 
is displayed above the diagonal, and genetic distance based as a percentage is displayed below the diagonal (A. betsileo was formerly 
referred to as Avahi species nova 41), 



1 

2 

"i 

4 

5 

6 

7 

8 

9 

1 


75 

68 

68 

74 

105 

238 

244 

247 

2 

3.21+0.36 


77 

31 

83 

(03 

225 

236 

242 

J 

2.78=0.36 

3.24*0.41 


69 

52 

107 

239 

241 

249 

4 

2.92*0.37 

1.33*0.23 

2.88*0.39 


74 

99 

220 

229 

235 

5 

3,14+0,37 

3.58*0.41 

2.09*0.32 

3.20*0.41 


110 

241 

246 

251 

6 

3.78*0,37 

3.66*0.35 

3.75*0.36 

3,45+0,37 

3.99*0.37 


246 

247 

251 

7 

11.06*0,68 

10.46*0.68 

11.06+0.69 

10.20*0.66 

11.29*0.71 

10,69+0.65 


67 

66 

8 

11.28*0,72 

10,96*0.71 

11.08=0.66 

10,59*0.68 

11.47*0.72 

10.69+0.69 

2.86*0.37 


64 

9 

11.44*0.69 

11.25*0,69 

11.48=0.63 

10.88*0.67 

11.72*0.68 

10,88*0.65 

2.83*0.32 

2.69*0.33 
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nova #1 to A. peyrierasi #3), respectively. The 
highest and lowest average uncorrected pairwise 
distances and absolute genetic differences for the 
PAST sequence alignment between the eastern 
woolly lemurs range from 3.99 to 1.33 and 110 to 
31 (A. m. meridionalis to A. laniger and A. species 
nova #1 to A. peyrierasi #3), respectively. The 
highest and lowest average uncorrected pairwise 
distances and absolute genetic differences for the 
D-Loop sequence alignment between the western 
woolly lemurs range from 4.17 to 2.76 and 24 to 
15 (A. unicolor to A. cleesei and A. unicolor to 
A. occidentalis), respectively. The highest and 
lowest average uncorrected pairwise distances 
and absolute genetic differences for the PAST 
sequence alignment between the western woolly 
lemurs range from 2.86 to 2.69 and 67 to 64 (A. 
unicolor to A. cleesei and A. unicolor to A. oc- 
cidentalis), respectively. 


The Bayesian analyses are presented in Figs. 
18-19. All analyses differentiate six recognized 
Avahi species, two types of A. peyrierasi and the 
newly described species. The minimum spanning 
network presents d^grammatically the speeiation 
among the six recognized woolly lemurs (Fig. 22). 
There is high bootstrap support for the ML analy¬ 
sis with respect to the topology among the genera 
and species (Figs. 18-19). The phylogenetic trees 
reconstructed through the various inferences and 
models produced identical topologies with the 
exception of the sister relationships between A. 
cleesei and A. occidental is or A, uni col or. Based 
on the sequence fragment utilized to reconstruct 
the tree; the D-Loop sequence fragment aligned A. 
cleesei with A. unicolor ? but the PAST sequence 
fragment aligned A, cleesei with A. occidentalis 
(D-Loop versus PAST fragments; Figs. 9-12). 


Discussion 


Notwithstanding what previously was consid¬ 
ered superficially inconspicuous phenotypic charac¬ 
ters, the woolly lemurs present subtle but definable 
morphological and phenotypic traits (Figs. 2-7,15-17, 
21-24; Table 1; Appendix Id). Milne-Edwards and 
Grandidier (1875a, b) noticed Avahi in the north¬ 
ern Ampasindava area in the late the 19th century 
and remarked that this Avahi from northwest was 
different and smaller than other vl laniger but did 
not recognize the significance of the differentiation 
at that time. The trend shown by the grouping of 
some morphological traits and the tendency towards 
lighter is evident in the data analyses (Tables 1A 
and IB). The tendency towards greater variation is 
expressed in the confidence intervals and the range 
of trait means observed more often among the eastern 
woolly lemur rather than in the western woolly lemurs 
demonstrated in Appendix Id. This should not be so 
surprising given the magnitude of variation in the 
mitochondrial haplotypes detected in the eastern as 
compared to the western species. Furthermore, the 
molecular sequence data submitted for the woolly 


lemurs corroborates the ^global” phylogeographic 
distribution of eastern and western lemur clades 
(Figs. 8-14; Louis et al. 2006a, 2006b; Zaramody et 
al. 2006). 

The D-Loop sequence presented was gener¬ 
ated for only a single representative for each recog¬ 
nized species of the family Indriidae, specifically 
the genus Propithecus (data were generated for two 
individual samples of Indri indri ; Figs. 8-10). The 
molecular sequence data offers a tentative glance at 
the relationships of this group to the other genera 
and species of lemurs. The relatively low bootstrap 
support (56%) between Indri indri and the genera 
Avahi and Propithecus should be further investi¬ 
gated, utilizing sequence data representing not only 
a coding region (for instance, PAST fragment or 
cytochrome B subunit gene), but also nuclear DNA 
sequence data and the addition of multiple samples 
from multiple subpopulations and/or localities 
that are representative of the recognized species’ 
distributions (Fig. 8). 


(text continued on page 43) 
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Figure 18. Bayesian majority rule maximum likelihood analysis cladogram derived from the D-loop DNA sequence 
data from the 77 Avahi individuals with 26 out-group taxa. Values above branches indicate support of bootstrap pseu- 
doreplicates with values greater than 50% reported. 
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Figure 19. Bayesian majority rule maximum likelihood analysis cladogram derived from the PAST DNA sequence 
data from the 77 Avahi individuals with 19 out-group taxa. Values above branches indicate support of bootstrap pseu- 
doreplicates with values greater than 50% reported. 
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Figure 21. Avahi meridionalis, southern woolly lemur, at Andohahela National Park (Manangotry). Photo 
by Edward E. Louis, Jr. 
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Figure 22. Avahi ramanantsoavana, Ramanantsoavana’s woolly lemur, at Manombo Special Reserve. Photo by 
Richard Randriamampionona. 
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Figure 23. Avahi betsilio, Betsilio woolly lemur, at Fandriana District (Bemosary Classified Forest). A. 
betsilio was formerly referred to as Avahi species nova # 1. Photo by John R. Zaonarivelo. 
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Figure 24. Avahi betsilio, Betsilio woolly lemur, at Fandriana District (Bemosary Classified Forest). A. 
betsilio was formerly referred to as Avahi species nova #1. Photo by John R. Zaonarivelo. 
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A comprehensive dataset that includes repre¬ 
sentatives from multiple populations and sequence 
fragments for all recognized genera and species will 
be required to conclusively resolve the taxonomic 
relationships between the family Indriidae and other 
genera, but this data does offer corroborative evidence 
for establishing species versus subspecies status for 
regional populations of the genus Avahi based on the 
phylogenetic species concept (Cracraft 1983). For in¬ 
stance, Zaramody et al. (2006) suggested that the status 
of A. cleesei was questionable based on a relatively low 
genetic distance between A. cleesei and A. occidenta- 
lis. By evaluating the uncorrected genetic distance, 
morphological, phenotypic, and pairwise aggregate 
analysis data, along with the inclusion of the third rec¬ 
ognized western Avahi , A. unicolor , we conclude that 
the parameters generated with the multiple data sets 
are analogous to those established with other species, 
including comparisons within the genera Propithecus 
and Lepilemur (Figs. 8-14; Tables 1-4; Appendix Id; 
Mayor et al. 2004; Thalmann and Geissmann 2005, 
2006; Andriantompohavana et al. 2006; Louis et al. 
2006a. 2006b). For example, there exist significant 
differences within the cytogenetic, phenotypic, and 
morphological data among sifaka species which con¬ 
sequently led to full species status for seven sifakas of 
the genus Propithecus (Mayor et al. 2004; Mittermeier 
et al. 2006 ). Despite the elevation to full species status 
for the sifakas, the genetic distance based on D-Loop 
sequence data between Propithecus tatterscilli and P. 
coquereli or P. edwardsi and P. diadema is relatively 
modest (1.79 and 1.96, respectively; Mayor et al. 2004). 
The minimal genetic distance between individual sifaka 
species, however, can establish supportive evidence in 
favor of the legitimacy of A. cleesei since the genetic 
distance between each pair of sifaka mentioned previ¬ 
ously is less than estimated between A. cleesei and A. 
occiclentalis (3.86) and/or A. unicolor (4.17; Mayor et 
al. 2004). 

Despite several citations which identify type 
locality of Avahi laniger in the Betanimena country, An- 
tongil Bay, Maroantsetra, the type locality of A. laniger 
remains unknown or uncertain (Jenkins 1987; Groves 
2001; Zaramody et al. 2006). Nonetheless, based on 
phenotypic and molecular data (PAST fragment only), 
two types of A. laniger from the extreme ends of its 
range are evident within these data (Figs. 1,4, and 15; 
Tables 2-4; Mantadia and Mananara-Nord National 


Parks). Further research, including comprehensive 
samples and data generation, will be needed, however, 
to conclusively resolve the taxonomic relationships of 
the woolly lemurs between the Mangoro/Onive Rivers 
and Anjanaharibe-Sud Special Reserve. 

The data also established two types of A. peyri- 
erasi (Figs. 8-14. 16-17; Tables 1-4). In forest frag¬ 
ments south and north of Ranomafana National Park 
(Tolongoina Classified Forest; Zaramody et al. 2006), 
unpublished data reveal that both types of A peynercisi 
are found within the same habitat (within Vohiparara, 
Ambatolahy Dimy, and Talatakely, all subpopulations 
within Ranomafana National Park). However, no be¬ 
havioral or hierarchal data exists within a population 
or family group, so further research (mitochondrial 
DNA sequence and genotype data (pedigree analysis) 
combined with long term behavioral observations) is 
required to establish the relationships between these 
two types of sympatric Avahi. Furthermore, we suggest 
that the two subspecies, A. meridional is meridional is 
and A. m. ramanantsoavana , described by Zaramody 
et al. (2006), should be elevated to full species because 
the molecular, phenotypic, and morphological data is 
equivalent, if not more conclusive, to the other previ¬ 
ously recognized woolly lemurs (Figs. 8-14, 18-22; 
Tables 1-4). 

In this paper, the current Avahi taxonomy was ex¬ 
amined according to the Phylogenetic Species Concept 
(PSC) sensu Wheeler and Platnick et al. (2000), Mayor 
et al. (2004), and Louis et al. (2006a, 2006b). The 
diagnostic characters or attributes define evolutionary 
significant units (ESUs). Several authors suggest that 
ESUs are equivalent to species as determined through 
the Phylogenetic Species Concept (Cracraft 1983; 
Barrowclough and Flesness f 996; Amato et al. 1998). 
The constant addition of samples to the PAA data set 
will continue to test the distinction of these characters. 
The identification of a new species in the following 
description establishes the essential need for extensive 
as well as detailed sample collections across Mada¬ 
gascar to determine geographic ranges and taxonomic 
status for all of the woolly lemurs. This research is 
especially needed on opposite aspects of all rivers (for 
example, north of the Mangoro River to Anjanaharibe- 
Sud Special Reserve). Due to the inability to maintain 
woolly lemurs as long-term live vouchers in captivity, 
whole blood, morphometric, and e-voucher photos will 
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serve as the type series for the newly described species, 
Avahi species nova #1, described below. In this case, 
an attempt was also made to identify existing museum 
specimens to represent the type series, but a museum 
specimen has not been identified from this specific 
region which could serve as the holotype. 

Avahi betsileo , New Species 

Type Series. —Whole blood for FAN2.14 (TK 
125757), adult male; FAN2.15 (TK 125758), adult 
female; and FAN2.20 (TK 125759), adult male; are 
stored and curated at the Museum of Texas Tech 
University. Individual measurements, e-voucher 
photos, and collection data are given in Avahi Field 
Data Appendix la. FAN2.14, FAN2.15, and FAN2.20 
were collected by Richard Randriamampionona, Jean 
C. Randriamanana, Gerard Nalanirina, Razafindraibe 
Jean, Rambinintsoa Andriantompohavana, and John 
R. Zaonariveloon 13 July 2002, 13 July 2002, and 15 
July 2002, respectively. 

Type Locality. —Madagascar, Province de 
Fianarantsoa, Region AmoronT Mania, District of 
Fandriana, Bemosary Classified Forest, 20°20'60.1 "S, 
47°3336.1 "E, and south of Mangoro River. 

Description .—When A. betsileo is compared 
to other southeastern woolly lemurs, Avahi betsileo 
(1.05 kg) is approximately the same size as A. pey¬ 
rierasi type #2 (0.86 kg) and A. meridionalis (1.06 
kg), but slightly larger than A. peyrierasi type #3 
(0.98 kg) and A. ramanantsoavana (0.98 kg; Table 
1). Phenotypically, A. betsileo differs significantly 
from other eastern woolly lemurs with a pelage that 
is primarily light reddish-brown on the body and on 
the dorsal surface of the extremities (Figs. 23-24). 
A. betsileo has a distinct facial mask with grayish 
pelage under the mandible and the ventral surface of 
the extremities. A. betsileo has diffuse cream colored 
eyebrow markings and a thicker pelage on the head 
which gives a more round or oval-like appearance of 
the head that distinguishes it from the other eastern 
woolly lemurs. The venter is dark gray towards the 
midline and diffuses to a light gray ventrolaterally. 
The tail is primarily reddish brown, darker on the 
dorsal surface than the ventral portion which is a 
lighter reddish blonde. 


Diagnosis. —In the D-loop and PAST sequence 
fragments, A. betsileo differs from the closest rela¬ 
tives, A. peyrierasi type 2, and A. peyrierasi type 3 by 
4.75%±0.91% (28 informative sites) - 2.92%±0.37% 
(68 informative sites) and 3.43%±0.77% (19 informa¬ 
tive sites)- 1.33%±0.23%(31 informative sites), respec¬ 
tively. In the D-loop and PAST sequence fragments, 
A. betsileo differs from the closest species relative to 
geographic distance, A. laniger , by 5.50%±0.91% (42 
informative sites) and 3.45%±0.37% (99 informative 
sites), respectively. 

Distribution. — A. betsileo is currentlv known 
in the Bemosary Classified forest (Fandriana) and the 
regional distribution will be tentatively set from south 
of the Mangoro River and north of the Mananjary River, 
until its distribution can be confirmed. 

Comparisons and Remarks.—Avahi betsileo 
is approximately the same size as A. peyrierasi type 
#2 and A. peyrierasi type #3, larger than A. raman¬ 
antsoavana, and smaller than A. laniger (Table l). 
Although A. betsileo is not genetically that different 
from both types of A. peyrierasi based on absolute and 
percentage differences, the phenotypic differences are 
noticeable (Figs. 23-24). Additional survey work is 
required to determine the southern and northern range 
of A. betsileo and the northern extent of A. peyrierasi 
type #2 and #3. 

Etymology. —The name betsileo is proposed for 
this species and is derived from the Malagasy tribe from 
the Fandriana region. 

Vernacular Name. —Betsileo woolly lemur. 

By applying the phylogeographic studies by 
Thalmann and Geissmann (2000) and Pastorini et al. 
(2003), geographic barriers can be defined that affect 
multiple taxa ultimately of both flora and fauna. By 
establishing and prioritizing these phylogeographic 
regions, conservation management decisions can be 
implemented that will support the preservation of bio¬ 
diversity within these defined regions. These areas of 
discrete or unique biodiversity 7 show that the all major 
rivers are significant factors or barriers and augment 
the processes of speciation (Fig. 25). The emerging 
data display a distribution of unique biodiversity 
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Figure 25. Distribution map of the woolly lemurs of Madagascar - Revised distribution of the woolly lemurs. The letter “A” 
represents the location of Ankarana Special Reserve. The letter “B” represents the location of Kalambatritra Special Reserve. 
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framing the topology between all rivers (Andrian- 
tompohavana et al. 2006; Louis et al. 2006a, 2006b; 
Rioux Paquette et al. 2006). This relationship be¬ 
tween the rivers and biodiversity transcends taxa 
and is quite remarkable and consistent (Louis et al. 
2006b). As discrete species from different genera 
are being described, their individual distributions 
can be overlaid, resulting in river bound “pockets” of 
unique biodiversity (Cover Figure; unique woolly and 
sportive lemurs forming paired units of distinct biodi¬ 
versity). The persistent and rapid loss of habitat and 
the resulting fragmentation of panmictic populations 
have compelled wildlife and conservation agencies to 
take protective action according to existing guidelines 
and information with the ultimate goal of prioritiz¬ 


ing species and/or sites. The explosive rate of the 
deforestation in Madagascar, however, has eliminated 
many of the available options (Green and Sussman 
1990). Because haplotypes can be unique to each 
population, simply preserving one population will not 
necessarily maintain species-wide genetic variability 
(Pope 1996; Louis et al. 2006b). Despite the recent 
revisions of the woolly lemurs, isolated fragments or 
forests that represent unique habitat in Madagascar, 
such as Kalambatritra and Ankarana Special Reserves 
and Masoala National Park, could potentially contain 
particular or distinctive Avahi species, so additional 
comprehensive field work and laboratory work needs 
to be completed (Cover Figure). 
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Appendix I 

The following Appendices to this publication are available online at the indicated website addresses. 

a. Avcihi Field Data Appendix (Individual data hie for each Avahi, including morphometries, photos, sequence 
accessions, global position system, microchip data, gender, and location). 

http://lQ.lQ.lQ.3/ccr/genetics/lemur/index.asp?page=ccr/genetics/lemur/Avahi.htm 

b. Appendix Primer Table I (Summary of designed primers for D-loop and PAST fragments). 

http://www.omahazoo.com/ccr/genetics/papers/appendixAvahiprimertablel.pdf 

c. Appendix Haplotype Table I (Summary of designated haplotypes for Avahi from Fandriana and 
Ranomafana localities for D-loop and PAST fragments). 

http://www.omahazoo.com/ccr/genetics/papers/AvahihaplotvpetableI.pdf 

d. Appendix ANOVA dendrogram of the morphological data summarized in Table 1. 

http://www.omahazoo.com/ccr/genetics/papers/appendixAvahimorphodata.pdf 


Sample (98 Avahi and 25 outgroups total) used in the present genetic study and taxonomic revision of the Madagascar lemur genus 
Avahi. a TK Number is referenced voucher curated at the Museum of Texas Tech University. b Mitochondrial ON A sequence data for 
D-Loop (D-Loop or control region) and PAST (Pasiorim fragment) for each sample are available from Gen Bank under the listed 
accession number. N/A signifies that sequence fragment is not available . (A. betsileo was formerly referred to as Avahi species nova 
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